On the Superradiant Phase in Field-Matter Interactions 

O. Castanos, E. Nahmad-Achar, R. Lopez-Pena, and J. G. Hirsch 
Instituto de Ciencias Nucleares, 
Universidad Nacional Autonoma de Mexico, 
Apdo. Postal 70-543 Mexico 04510 D.F. 

Abstract 

We show that semi-classical states adapted to the symmetry of the Hamiltonian are an excellent 
approximation to the exact quantum solution of the ground and first excited states of the Dicke 
model. Their overlap to the exact quantum states is very close to 1 except in a close vicinity of the 
quantum phase transition. Furthermore, they have analytic forms in terms of the model parameters 
and allow us to calculate analytically the expectation values of field and matter observables. Some 
of these differ considerably from results obtained via the standard coherent states, and by means 
of Holstein-Primakoff series expansion of the Dicke Hamiltonian. Comparison with exact solutions 
obtained numerically support our results. In particular, it is shown that the expectation values 
of the number of photons and of the number of excited atoms have no singularities at the phase 
transition. We comment on why other authors have previously found otherwise. 



PACS numbers: 42.50.Ct, 03.65.Fd, 64.70.Tg 



I. INTRODUCTION 



The field- matter interaction is commonly described by the Dicke Hamiltonian [TJ, which 
considers N two-level atoms with energy separation equal to Hua immersed in an external 
one-mode electromagnetic field of frequency Qp inside a cavity. This Hamiltonian comes from 
the standard quantization of the multipolar form of many atoms interacting with classical 
radiation [2J in the long wave approximation. The Dicke model was established exploiting the 
analogy with the spin. It is considered a multiatom generalization of the Jaynes-Cummings 
model [3J, a completely soluble quantum Hamiltonian of a two-level atom in a one mode 
electromagnetic field. The Tavis-Cummings model [I] is also a many body extension of the 
Jaynes-Cummings model, which considers the rotating-wave approximation. 

A very important feature of the Dicke Hamiltonian discovered by Hepp and Lieb [5] is the 
presence of a phase transition from the normal to the superradiant behaviour. The latter, 
introduced by Dicke to describe coherent radiation, is a collective effect involving all N 
atoms in the sample, where the decay rate is proportional to iV 2 instead of N, the expected 
result for independent atom emission. While this transition to the superradiant phase has 
been much debated in the literature [EH9], due to the smallness of the field-matter coupling 
strength, recent experimental results indicate that it can actually be observed [T014T2] . 

In this contribution we use the tensorial product of SU(2)- and if iy(l)-coherent states 
as a trial state to calculate the expectation value of the Hamiltonian and construct its 
energy surface. This is a function that depends on four variables (q,p,8,<f)) which define 
the coherent states, and two essential parameters (um,7) which are associated to the energy 
separation of the two-level atoms and the coupling strength between the atoms and the field, 
respectively. By studying the stability properties of the critical points of this energy surface 
we obtain the best variational wavefunction and a separatrix which divides the parameter 
space in two regions: normal and superradiant behaviour, according to 7 < 7 C or 7 > j c , 
respectively, with 7 C = ^JujaJA. 

The Dicke Hamiltonian is invariant under transformations of the cyclic group C2, a sym- 
metry which is not preserved by the trial state. The symmetry can be restored by acting 
with the projectors of the symmetric and antisymmetric representations of the cyclic group 
C*2- We determine the statistical properties of field and matter observables in the restored 
variational states \^p) in the superradiant region. The comparison with the same observ- 



2 



ables evaluated from the exact diagonalization of the Hamiltonian at finite number of atoms 
TV exhibits an excellent agreement between them, which improves as N becomes larger [T3] . 
The fidelity between both types of states is also evaluated; we establish that the fidelity in 
the superradiant regime is very close to unity, dropping fast in a small neighbourhood of the 
separatrix. The great advantage of working with the restored-symmetry states is that they 
allow us to write down analytic expressions for all the expectation values of field and matter 
observables, as well as with important mathematical relations. These analytic expressions 
are exact in the thermodynamic limit, when the number of particles goes to infinity. 

The Dicke Hamiltonian has been studied thoroughly by Emary and Brandes [H] with 
a variational procedure in terms of two if iy( Incoherent states. One of them is related to 
the radiation field degrees of freedom and the other arises from a Holstein-Primakoff (H- 
P) realization of the SU(2) generators [15]. They obtain the appropriate relation between 
the parameters of the Dicke Hamiltonian to describe the phase transition, as well as a 
description for the normal phase in the thermodynamic limit. However, they consider in the 
superradiant region, an approximated Dicke Hamiltonian coming from a series expansion of 
the H-P realization truncated to second order in terms of the ratio between the number of 
excited atoms over the total number of atoms, assuming that is a very small quantity. Under 
this assumption Nagy et al. [11] found singularities at 7 = 7 C in the expectation values of the 
number of photons and the number of excited atoms, and related them with the divergences 
in their corresponding fluctuations [H] . 

In this contribution we extend and provide details of the calculations presented in [13] . 
The expectation values of the mentioned observables are calculated by means of the symme- 
try adapted variational states and find that there is in fact a divergence of the expectation 
values, but that this remains throughout the superradiant region, contrary to the published 
results: there is no singularity. The expectation value of the number of photons (a)a) and 
the squared fluctuation of the first quadrature of the electromagnetic field (Ag) 2 are propor- 
tional to the number of atoms at the phase transition and in the superradiant regime, while 
the squared fluctuation (AJ X ) 2 is proportional to N 2 in the same region. Comparison with 
the exact numerical solution supports our results. This could have consequences in studies 
of entanglement and chaos in the superradiant region. 
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II. DICKE HAMILTONIAN AND SYMMETRY ADAPTED STATES 



The Dicke Hamiltonian involves the collective interaction of N two-level atoms with 
energy separation HO a with a one- mode radiation field of frequency up in the long wavelenght 
limit. It has the form 

H D = a) a + u A J z + -y= (a f + a) ( J+ + , (1) 



'N 

where oja = &a/&f is given in units of the frequency of the field, and 7 = ^/Cj f = 
a/ {2-ng/hCop) dba-ep, is the coupling parameter between the matter and field. The g denotes 
the density of atoms in the quantization volume, d\, a is the excitation matrix element of the 
electric dipole operator of a single atom, and e*p the polarization vector. The operators a', a 
denote the one-mode creation and annihilation photon operators; J z the atomic relative 
population operator; and J± the atomic transition operators. We will consider completely 
symmetric states for which we have N = 2j. 

Later we will find it convenient to divide the Dicke Hamiltonian by the total number 
of particles, having in this way an intensive Hamiltonian operator. In order to study the 
thermodynamic limit one would take N — » 00. 

To find the symmetries of the Hamiltonian, one considers the unitary transformation 



U = exp ^200^ with A = a) a + J z + yi 2 + 1/4 - 1/2 denoting the excitation number 
operator. One can show that 

(J a C/t = e-^° a , UJ + U j = e"^ J+ , (2) 

with the corresponding hermitean conjugated relations. Substituting these results into the 
expression for the Dicke Hamiltonian one has 

UH D U j = d ] d + u A J z + -j= (a ] L + a J+) + (V 2 ^ a f J + + e 2i(t>0 a J_) . (3) 



/N V / VN 

If the counter-rotating term (last term of the expression) is neglected, one recovers the Tavis- 

Cummings model and it is immediate that it is invariant under any rotation with arbitrary 

0o- Their symmetrized solutions have been thoroughly studied in [161 Ej- F° r the Dicke 

Hamiltonian, it is necessary to restrict to rotations by an angle 0o = 0, tc for the Hamiltonian 

to remain invariant. Thus the invariance group for the Dicke model is C2 

The projection operators for this group are 



I±e twA ) , (4) 
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which allow us to write the eigenf unctions of the Hamiltonian displaying its symmetry 
explicitly. In terms of the eigenstates \v) of the photon number operator, and the square of 
the angular momentum and its projection in the z-axis \j, m), they may be expressed as 

oo A 

I4t> = £ E (5) 

A=0 i/=max{0,A-2j} 

where contains only even values of A = v + j + m, the eigenvalues of A, while \</>j.-) 

contains the odd ones. The index k denotes the state number of the even and odd solutions. 
In practice one uses a maximum value for A which guaranteed the convergence in the energy 
eigenvalues of the ground and first excited states. The coefficients c k x v can be obtained from 
the diagonalization of the Hamiltonian matrix, whose dimension is 

^ J 2 ("^ max "0 (^max 2) , A max ^ 2j , 

[ (2j + 1) (A max - j + 1) , A max > 2j , 

in the basis states without definite parity of A. This matrix, in the symmetry adapted basis 
states, breaks into two pieces of even and odd parity respectively. For A > 2j, the dimensions 
of the matrices are 

d + = (j + l) 2 + (2j + l)(s + - j) , = j(j + 1) + (2j + - j) , 

for integer j, with s + = \\ ma x/2\ and s_ = [(\ m ax + 1)/2J . For A < 2j, they are d+ = 
(s + + l) 2 and d- = + 1). Similar expressions can be obtained for half integer j. 



III. ENERGY SURFACE AND CRITICAL POINTS 

The energy surface is found by taking the expectation value of the Dicke Hamiltonian with 
respect to the tensorial product of coherent states of Heisenberg-Weyl and SU{2) groups 
|a) ® 10 , given by 



oo 

a 

\at 



) = exp(-|a| 2 /2) Y.^=>) 
10 = 1 ' y ( 23 V^C^Ii, rn) 



m=-j 

where the parameters a and ( are complex numbers. Using the Ritz variational principle 
one may find the best variational approximation to the ground state energy of the system 
and its corresponding eigenstate. 



To calculate the expectation values of the field and matter observables it is convenient 
to find the representation of the angular momentum and Weyl generators with respect to 
the tensorial product |a) <g> |C) [IB]- For example, for the annihilation photon operator we 



consider the matrix element 



M 2 

(a\a\ip) = e" 2 (0|e a a a 
d a 



+ o <<#> > ( 6 ) 



da* 2 

where \ip) denotes an arbitrary state. We proceed in a similar form for the other observables. 
The representations of the annihilation and creation photon operators are thus given by 

d oc ^ + * 
fl ^ + 2' a ^ a ' 

and for the matter observables we have 

j _> i— + C — 

z 1 + lCl 2 d(*' 

I Y * 2 + ^l 2 +(C*) 2 — 

1 + ICI 2 + u j «9C 

" 1 + ici 2 + oc* ■ 

The energy surface can be then calculated straightforwardly (TUl E] 

H(a,() = (a\®((\H D \a)®\C) 

= -(p 2 + q 2 ) — j uja cos 6 + 2yj7 g sin 6 cos . (7) 

In this expression we take the harmonic oscillator realization for the field part and the 
stereographic projection for the angular momentum part, 

1 6 

a = —(q + ip) , C = e _i *tan-, (8) 

where (q,p) correspond to the expectation values of the quadratures of the field, and (#,</>) 
determine a point on the Bloch sphere. 

The minima and degenerate critical points are obtained by means of the catastrophe 
formalism [2U]. By calculating the Hessian of the energy surface, we see that when 7^ = Cou/4 
the critical points degenerate and, for that value of the field-matter coupling, the phase 
transition from the normal to the superradiant behaviour of the atoms takes place. 
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The critical points which minimize H are given by 

6c = 0, q c = 0, p c = 0, for |7| < 7c, 

6» c = arccos(7 c /7) 2 , q c = -2 ^ j y/l - (7^7) 4 cos C , p c = , for (7! > 7 C . 

The last column of this array shows the conditions in parameter space to guarantee that they 
constitute a minimum of the energy surface. The first row describes the minimum critical 
points for the normal phase, while the second row describes those for the superradiant 
regime. In this last case one has <j) c = 0, 7r. (If one were to consider a case in which ua < 0, 
we would then have 7^ = — uja/4 and the critical points for the normal phase would be 
9 C = it, q c = 0, p c = 0, for |7| < 7 C ). It is convenient to work with the variable x = 7/7 c in 
terms of which the energy values for the minima just described are 

-^normal ~~ — 2A^7 C , -^superradiant ~~ — ^ 1c x (l + x ) • (10) 
In a similar form, the expectation values of A at the minima are 

^normal — ^ > ^superradiant — T — x + ^7 C a; (1 — x )] . (11) 

The fluctuations AA are zero in the normal region, and in the superradiant phase take the 
form 

AA = ^/f Q + 2 7 c 2 ^)(l-x-4) . (12) 

To study the statistical properties of the variational states in the superradiant regime we 
calculate the expectation values of linear matter and field observables with respect to the 
tensorial product of coherent states, as well as their fluctuations. The results are given in 
the left side of Table [Tj Notice that we have introduced the quadratures operators q and p of 
the electromagnetic field and they can be written in terms of the creation and annihilation 
photon operators. 

IV. SYMMETRY-ADAPTED COHERENT STATES 

To build variational states which preserve the symmetry of the Dicke Hamiltonian we 
apply the projectors of even and odd parity |4]) to the coherent states |a) <8> |C), to obtain 

k C>± = A4 (\a) ® |C) ± I - a) ® I - 0) , (13) 
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where the normalization factors M± are given by 



A4- 2 = 2(l±exp(-2H 2 )('i T |P) JV ) . (14) 



The expression ( 13 ) is useful because it allows us to calculate the expectation values of the 
relevant operators in closed form, in particular the expectation value of the Hamiltonian. 
Evaluating the energy surface with respect to the symmetry adapted states, one obtains 



{H)± = 4^ + ^{ 1 -l±e^)(cos^} 

N { , , tan 2 # cos 9 

ua s (cos 9) ± 



2 ' l±e±(p 2 +9 2 )(cos^) 



/— — I ±p tan 9 sin 6 + qe p2+q2 sin 6* cos0 (cost , , 
+ ^ 7 { e^(oo S 6>)-"±l j ■ (15) 

In the limit N — >■ oo, (H) + reduces to Eq.([7]), when | cos^l ^ 1. However, it is important to 
stress that the analysis can be carried out for any value of N, and in this contribution we 
work at finite N. 

The traditional approach in many-body physics is to use the critical points in Eq.Q of 
the original energy surface, in order to obtain the trial state which approximates the two 
lowest energy states, and in which to evaluate the expectation values of the observables. 



-N 



Formally, one should calculate the critical points of (15), instead of that, we use the critical 
points associated to Eq.([7]). However all of their critical points also are extreme values for 
(15), except in a small vicinity of the separatrix, as shown in Fig.Q. In that figure, the 



derivatives of the energy surface, given in (15), with respect to q and 9, evaluated at the 



critical points q c and 9 C , are plotted against 7. (The derivatives with respect to the other 
two variables, p and 0, are identically zero). The neighborhood around the separatrix, where 
they differ, diminishes as j increases (cf. Fig.Q). For the superradiant region in which the 



points are critical , they are also minima. The stability analysis of the energy surface (15) 
and the study of their behaviour in the near-neighborhood of the separatrix will be reported 
elsewhere 



By substituting the critical points p c = 0, q c = —2 y/j ^ c x \Jl — x 4 cos0 c , and cos# c 



x in (|15|), the energy surface associated to the superradiant regime takes the form 



(H) ± = -iV 7c V 



x 



1±F 



(16) 
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where we have defined 



x 



F = x 



-2N -2N 1 lx 2 {l-x- i ) 



(17) 



The overlap between the ordinary coherent states and our symmetry-adapted states can be 
written as 

(18) 



d{H) 
dq c 



0=io) 



■ 
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FIG. 1: Derivatives of the energy surface of the symmetry-adapted states, evaluated at the critical 
points q c and 9 C , are plotted as a function of the interaction strength 7. We use the frequency 
loa = 1 and N = 20 atoms. Continuous lines refer to the ground state while dashed lines to the 
first excited state. 

In Fig. Q the energy surfaces for the even and odd parities of A are shown together 
with those for the exact solution, showing a remarkable agreement between them in spite 
of having considered a small number of atoms iV = 20. For the odd parity case, we have 
proposed, in the normal region (I7I < 7 C ), a combination of states with A = 1: 



)i> = cosfi |0) ® \j, -j + 1) - 3- sinfi |1) <g> \j, -j) 

7 



(19) 



and minimized the expectation value of the Hamiltonian with respect to Q. One finds 
tan(2f2 c ) = 2 |7|/(1 — w^). This gives, for the resonant case, Q c = it /A. 

V. STATISTICAL PROPERTIES 



A. Expectation values 



One can calculate the expectation values of the matter and field main operators in the 
superradiant regime with respect to the symmetry adapted states. The results are collected 
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dq 




FIG. 2: Derivative of the energy surface of the adapted-symmetry states, calculated at q c and 9 C , 
is displayed as a function of the interaction strength 7, for different number of atoms. The plots 
correspond to N = 20, 50, 100 atoms. It can be seen that the neighborhood of discrepancy from 
zero diminishes as N increases. In all cases ua = 1- 
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FIG. 3: Energies of the ground and first excited states of the Dicke Hamiltonian as a function of 
the interaction strength 7, for the frequency uja = 1 and N = 20 atoms. The thin and thick con- 
tinuous lines show the ground and first excited energies, respectively, evaluated with the projected 
variational states, while the dashed and dotted lines correspond to the ground and first excited 
energies calculated through the diagonalization of the matrix Hamiltonian. 

in Table [T] together with those associated to the coherent states. Both are evaluated at the 
critical points for the energy surface, and are given in terms of x and F. 
Observe that, for expectation values different from zero in the symmetry adapted states, 
the coherent state results can be obtained from the former by letting F go to zero, with 
the exception of (Ag) 2 and (AJ X ) 2 . These two exceptions arise because (q) = (J x ) = for 
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TABLE I: Expectation values and fluctuations of matter and field observables for the coherent 
and symmetry-adapted states in the superradiant regime. The mean-field behaviour obtained in 
the normal region can be recovered by taking the limit x — >■ 1. 





Coherent 


Symmetry Adapted 
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i T 2iV 7 2 x 2 (^)F 
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' 1± (^v-i)(i-^ 4 )^\ 


(AJ Z ) 2 


f (1 " ^ 4 ) 


T [lT(iV-l)(l-x 4 )F-x 4 


F 2 ] 


(Aata) 2 




iV 7 2 x 2 (I-*" 4 ) 


iV 7c 2 x 2 |iV 7c 2 x-6 (l-x 4 )( 7 i|) 2 










+ (1 - x- 4 ) 


"iV 7 2 x 2 (i_x- 4 ) + 7 i|" 


} 


{J z a)a) 


-iV 7c 2 (I-*" 4 ) 


-iV7 2 


x 4 (l-x" 4 ) (^) 




0x0) 


-^7c*(l-*- 4 ) 




V 2 ' cX 1±F 





the former but not for the latter. That F tends very quickly to zero as a function of 7 (or 
x) can be seen in Fig. (4), especially for large N. This is why coherent states have been 
so successful in the past as trial functions. It can be seen that the quantities which differ 
the most are some fluctuations of matter and field observables. To assess this difference we 
evaluated the fluctuation of the atomic transition operator J x and of the first quadrature of 
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the electromagnetic field, q. 




1.2 1.4 1.6 1.8 



FIG. 4: The behavior of F as a function of x is presented for a number of atoms in the cavity 
equal to N = 2, 10, 20 and 100. We use 7 C = 0.5 (in resonance). 



B. Fluctuations AJ X and Ag 

In Fig. (JsJ) we show the fluctuation of the transition operator (AJ X ) 2 for the ground and 
first excited states of the Dicke model with N = 10 atoms and oja = 1, calculated for the 
symmetry adapted states and the coherent states. While (AJ X ) 2 /N 2 tends erroneously to 
zero for the coherent state, as x increases, the results for the symmetry-adapted states show 
the appropriate quadratic dependence in N, asymptotically reaching (AJ X ) 2 w N 2 /4. In 
Fig. ^ we show (Aq) 2 for the ground and first excited states of the Dicke model with 
N = 10 atoms and ua = 1 calculated for the symmetry-adapted states and the coherent 
states. While for the coherent state the result is a constant 1/2, for the symmetry adapted 
states one has a linear dependence in the number of atoms. Additionally, in both cases one 
can see that the symmetry-adapted states compare really well with the exact result obtained 
from the diagonalization of Hamiltonian for the even and odd parity states, except in the 
close vicinity of the separatrix of the physical system [13]. These results show clearly the 
benefit of using the symmetry-adapted states over the coherent ones. 

C. Joint probability distribution function 

The joint probability of finding v photons and n e = j + m excited atoms, for even and 
odd parity states, is obtained by taking the modulus square of the scalar product of the 
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FIG. 5: Fluctuation of the atomic transition operator J x evaluated for the ground and first excited 
states of the Dicke Hamiltonian as a function of the interaction strength 7, for the frequency uja = 1 
and N = 10 atoms. The dashed dark and light lines show the fluctuations for the ground and first 
excited states, respectively, evaluated with the projected variational states, while the continuous 
lines describe the correspodning fluctuations calculated through the diagonalization of the matrix 
Hamiltonian. The result for the coherent state is shown in the dashed-dotted line. (Color online.) 



Fock and angular momentum states \v) <g> \N, n e ) with the symmetry-adapted states. The 
result depends on the values of N, j c , and x, and it is given by 

1 



V±(v, n e 



[1±("1H ^[X-t* ( l ~^)Y 

-2 \ n e 



X 



1 



X 



1 + x- 



2\ N-n e 



X 



N 



1±F 



(20) 



The joint probability distribution functions, for the even- and odd-parity states, are shown 
in Figs. ([7j [8]) for a system of iV = 10 atoms and in resonance, i.e., ua = 1. In Fig. ([7]), we 
consider 7 = 0.55 and in Fig. (|8p, 7 = 1.0. It is clear that there are holes in the distributions 
of each plot, corresponding to the forbidden values (odd or even) of the sum v + n e . Note 
that the Poissonian distribution when near the separatrix (7 = 0.55) turns into a quasi- 
normal distribution as we move away (7 = 1). While very few photons and excited atoms 
contribute to the state in the first case, many more contribute in the second case, as is to 
be expected. 

By summing over the number of photons v or the excited atoms n e , one determines the 
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FIG. 6: Fluctuation of the electromagnetic field first quadrature operator q evaluated for the 
ground and first excited states of the Dicke Hamiltonian as a function of the interaction strength 7, 
for the frequency oja = 1 and N = 10 atoms. The dashed dark and light lines show the fluctuations 
of the ground and first excited states, respectively, evaluated with the projected variational states, 
while the continuous lines display the fluctuations of the ground and first excited states calculated 
through the diagonalization of the matrix Hamiltonian. The result for the coherent state is shown 
in the dashed-dotted line. (Color online.) 



corresponding marginal distributions. For the distribution of the number of photons one has 

(21) 



±F 



while for the distribution of the number of excited atoms the result is 



V±(n e 



x~ 



2\ n e 



-1 , — 9\ N—n e -, 
1 + X Z \ 1 ± 



-ir x 2n f 



1±F 



(22) 



The behavior of these probability distribution functions for even and odd parity states is 
shown in Figs. ^ and (10) for values of 7 close to (x = 1.1) and far from (x = 2.0 ) the 
separatrix. For x = 1.1, they are different, while for x = 2, they cannot be distinguished. 
We notice that the binomial and Poissonian distributions become quasi-normal as we move 
away from the separatrix. This is to be expected, since, from the results given in Table it 



is immediate that the distribution function of photons, Eq. (21), can be rewritten as 

'1 ± (-lYx~ 2N ' 



V+(v) 



1 



1 ± e- 2 ^x~ 2N 



(23) 
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FIG. 7: The joint probability distribution function of photons and excited atoms is shown, for 
7 C = 0.5, 7 = 0.55, and N = 10 atoms. At the left, we plot the result for the even parity states 
while at the right we indicate the distribution for the odd ones. 



vn e + r vn e - 




FIG. 8: The joint probability distribution function of photons and excited atoms is shown, for 
7 C = 0.5, 7 = 1.0, and N = 10 atoms. At the left, we plot the result for the even parity states 
while at the right we indicate the distribution for the odd ones. 

where /x = (a)a) is the expectation value of the photon number operator in the coherent 
state. For large values of N this distribution becomes a Poisson distribution, which in this 
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limit is equivalent to the normal distribution 

1 



P±(u) 



cxp 



2[L 



(24) 



In Fig. ([9]), the distribution of the number of excited atoms for the even- and odd-parity 
states are shown. For x = 1.1, close to the separatrix, they are different, while for x = 2, 
we can not distinguish them. 
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FIG. 9: The probability distribution of excited atoms is shown, for 7 C = 0.5 and N = 10 atoms. 
At left we plot the result for 7 = 0.55, while at the right for 7 = 1. In both cases the red bar (left) 
indicates the distribution of the even parity states and the green bar (right) the result for the odd 
ones. 
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FIG. 10: The probability distribution of photons is shown, for j c = 0.5 and N = 10 atoms. At 
left we plot the result for 7 = 0.55, while at the right for 7 = 1. In both cases the red bar (left) 
indicates the distribution of the even parity states and the green bar (right) the result for the odd 
ones. 
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For the distribution function of the excited number of atoms, in the limit when iV 3> 1 
one has 

(//, - ( i - ./• -) .\ n' 1 

exp 



V±{n e 



(1 - x- 4 ) N/2 



(l-x~ i )N/2 



which for large values of x also gives a normal distribution: 



V±{n e 



1 



7T 



N/2 



exp 



K - N/2)' 
N/2 



(25) 



(26) 



VI. PROJECTED STATES VS. EXACT SOLUTION 



By substituting the critical points of the superradiant phase into the symmetry adapted 



variational states (|13| 14|), we obtain 

exp [- 



l^sr; Csr): 



±1 



^/2 N + l (1 ± F) 

1/2 

(1-x- 2 ) 



N 2 2 



. N 7 C x 

(i-^i^tiiHr)^- 



2 \ (v+n e )/2 



(l + x- 2 ) 



2 \(N+u-n e )/2 



v) ® 



N 
2 ' 



n, 



>■ 



(27) 



In the limit i > 1 we have F — > 0, a sr — >■ v#7 c i, and Csr — ^ 1- This implies that each 
matter coherent state goes to an eigenfunction of J x , i.e., a rotation by | acting on the state 
|j, — j). In that limit the Dicke model becomes completely integrable [13] . 

A good measure of the distance between quantum mechanical states is given by the 
fidelity; for pure quantum states it measures their distinguishability in the sense of statistical 
distance [22], but it is customary to use the fidelity as a transition probability regardless of 
whether the states are pure or not. 

In this contribution we calculate the fidelity of the exact ground and first excited states 
with respect to the corresponding symmetry adapted states 



I (V>P Inexact) I 2 • 



(28) 



In both cases the fidelity gives a result very close to 1, except in the vicinity of the quantum 



phase transition. In Fig. (11 ), the fidelity is shown as a function of 7 for iV = 10, 20, 40 and 



50. For increasing iV the fidelity rises more sharply to 1 as we move away from the phase 



transition, located at 7 C = 0.5. Note that (cf. Eq.(18)), had we used the coherent states as 
trial functions, we would at best obtain a fidelity F close to 1/2, in the limit of large 7. 
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0.7 0.8 0.9 1.0 




FIG. 11: Fidelity between the symmetry- adapted and exact quantum states, as a function of the 
interaction strength 7, shown for the ground state (even parity, left) and the first excited state 
(odd parity, right). In both cases, the plots show that for increasing N = 10, 20, 40 and 50 the 
fidelity rises more sharply to 1 as we move away from the phase transition (at j c = 0.5). In all 
cases da = 1. (Color online.) 

VII. DISCUSSION AND CONCLUSIONS 



By means of the symmetry properties of the Dicke Hamiltonian, we found analytic ex- 
pressions for the ground and first excited states, which allows also to determine in closed 
form the expectation values of matter and field observables. The procedure was the follow- 
ing: we calculated the expectation value of the Hamiltonian with respect to the tensorial 
product of Weyl and SU(2) coherent states. This expectation value defines a function called 
energy surface depending on phase space variables and parameters. We then determine the 
degenerate and minimum critical points of the energy surface. The minimum critical points 



yield the expressions (10) which give the mimimum energy of the system together with 



information about the quantum phase transition present in the Dicke Hamiltonian. This 
quantum phase transition occurs when u^up = A^ 2 . Afterwards we restore the symmetry 
exhibited for the Dicke Hamiltonian by means of the projection of the states to definite 
parity of the excitation number operator. These symmetry considerations lead to determine 
the ground and first excited states together with important differences from the mean field 
results for some fluctuations of the matter and field observables. Additionally they provide 
us with the joint probability distribution functions of the number of photons and excited 
atoms. 

The condition to transit from the normal to the superradiant regimes is very difficult 
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to satisfy for optical systems because the available dipole coupling strengths are usually 
smaller than the transition frequency. Some proposals to overcome these problems have been 
discussed in [TT], where it is stated that the quantum motion of a Bose-Einstein condensate 
trapped in an optical cavity can be used to realize the Dicke model. Indeed, this has been 
reported in [12J, where the Dicke Hamiltonian has been physically realized in a superffuid 
gas moving in an optical cavity. However, the feasibility of reaching the superradiant phase 
of the field- matter interaction is still under debate [23], |2l] . 

The predicted A/iV-behaviour for the collective iV-atom interaction strength of the model 
has been observed experimentally in circuit QED for up to 5 qubits [25J. 

In this work we have shown that semi-classical states adapted to the symmetry of the 
Hamiltonian are an excellent approximation to the exact quantum solution of the ground and 
first excited states of the Dicke model in the superradiant phase. Their overlap to the exact 
quantum states is very close to 1 except in a close vicinity of the quantum phase transition 



(cf. Fig. (11)), whereas that of the ordinary coherent states would be at best equal to 1/2 



(cf. Eq.(18)). Our projected states have analytical forms in terms of the model parameters 
and allow us to calculate analytically the expectation values of field and matter observables. 
We have found that in the superradiant regime the fluctuation (Ag) 2 of the first quadrature 
of the electromagnetic field is different from the results obtained via the standard coherent 
states and its value grows as a linear function of N. Something similar happens for the 
fluctuation in the dipole transition operator (AJ X ) 2 , where one finds a quadratic dependence 
in the number of atoms. Both these results contradict those obtained previously [H]. The 
expectation values of the number of photons and of the number of excited atoms were 
also studied, finding that there are no singularities at the phase transition. Those found 
previously [11] are an artifact of an inappropriate truncation of the Hamiltonian. The joint 
probability distribution functions for the ground and first excited states were shown, which 
may be used to characterize the states of atoms in a cavity according to whether n e + v has 
an even or odd value. 

In a future contribution, we will evaluate other properties like the entanglement entropy 
between field and matter, and the squeezing parameter for the electromagnetic field and 
atomic components, which has been the subject of much interest citevidal. 

Finally we want to remark that the present formalism allows for a simplification of the 
exact quantum calculation by considering an expansion of Dicke and photon number states 
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running from a minimum value \ m in to X ma x, where these values can be estimated from A c 
and its corresponding fluctuations 5X C as functions of 7, the coupling parameter between 
the field and matter. The same simplification can immediately be done in the proposed 



state ( 13 ) without changing the obtained results for all the expectation values and probability 
distributions. 
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